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ABSTRACT 

The cosmic star formation rate is observed to drop sharply after redshift z = 2. We 
use a large, cosmological, smoothed particle hydrodynamics simulation to investigate 
how this decline is related to the evolution of gas accretion and to outflows driven 
by active galactic nuclei (AGN). We find that the drop in the star formation rate 
follows a corresponding decline in the global cold-mode accretion rate density onto 
haloes, but with a delay of order the gas consumption time scale in the interstellar 
medium. Here we define cold-mode (hot-mode) accretion as gas that is accreted and 
whose temperature has never exceeded (did exceed) 10^'^ K. In contrast to cold- 
mode accretion, which peaks at z « 3, the hot mode continues to increase to z k, 1 
and remains roughly constant thereafter. By the present time, the hot mode strongly 
dominates the global accretion rate onto haloes. Star formation does not track hot- 
mode halo accretion because most of the hot halo gas never accretes onto galaxies. 
AGN feedback plays a crucial role by preferentially preventing gas that entered haloes 
in the hot mode from accreting onto their central galaxies. Consequently, in the absence 
of AGN feedback, gas accreted in the hot mode would become the dominant source 
of fuel for star formation and the drop off in the cosmic star formation rate would be 
much less steep. 

Key words: stars: formation - galaxies: evolution - galaxies: formation - intergalactic 
medium - cosmology: theory 



1 INTRODUCTION 

The cosmic star formation rate (SFR) is observed to peak 
between redshifts 2: « 3 and z 2 after which it decrease s 
by an order of magnitude (e.g. iHopkins fc BeacomI l20o3 ). 
Its evolution is thought to be determined by the combina- 
tion of the formation and growth of (dark matter) haloes, 
which depends on cosmology, and the distribution of SFRs 
in haloes as a function of halo mass and redshift. The latter 
depends on processes such as gas accretion, stellar mass loss, 
radiative cooling, (re-)ionization, and feedback from star for- 
matio n and active galactic nuclei (AGN) (e.g. ISchave et al.l 
I2OIOI ). 

The observed rates of star formation in galaxies can 
only be sustained for long peri ods of time if the galaxies 
are being fed continuously (e.g. iBauermeister et"al.i.2010.) . 
This feeding happens through the accretion of gas from the 
intergalactic medium. In the simplest picture of spherical 



E-mail: fvdvoortiOlstrw.Ieidenuniv.nl 



collapse, it is assumed that all gas accreting onto a halo is 
shock-heated to the virial temperature of the halo, reaching 
a quasi-static equilibrium supported by the pressure of the 
hot gas. If the cooling time of the halo gas is sufficiently 
short, then it may subsequently enter the central galaxy. 
Whether or not the halo gas can accrete onto the galaxy 
therefore d epends on both the tem perature and the density 
of the gas (|Rees fc Ostrikeilll977^ . However, if the cooling 
time of gas that has gone through a (hypothetical) accretion 
shock at the virial radius is short compared with the age 
of the Universe, then there can be no hydrostatic halo and 
hence a lso no virial shock. Th i s will be th e case for low-mass 
haloes (|Rees fc Ostrikeij|l977l : IWhite fc Reeslll978i ). The ac- 
cretion rate onto the central galaxy t hen depends on the in - 
fall rat e, but not on the cooling r ate l|White fc Frenklll99ll ). 
Indeed, iBirnb 01m fc Dekell (|2003D showed that a virial shock 
is absent for low-mass haloes in a spherically symmetric sim- 
ulation. 

However, cosmological simulations show significant de- 
viations from spherical symmetry. Galaxies form inside the 
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filaments that make up the cosmic web and these filaments 
contribute signi ficantly to, or eve n dominate, the gas sup- 
ply of galaxies l|Dekel et al.l 120091 '). Because the density in- 
side filaments is higher than that of the ambient gas, the 
cooling time of the gas is much shorter and the gas can 
accrete cold onto haloes of higher mas ses than predicted 
by spherically symmetr ic models (e.g. iKeres et al.l I2OO5I : 
iDekel fc BirnboimllioO^ ). This cold gas can reach the high 
densities of the interstellar medium (ISM) much more effi- 
ciently than the tenuous, hot gas in the halo. It is this 'cold 
mode' of accretion tha t predominantly feeds the galaxy and 
powers star formation feeres et al ]|2009l :[ Brooks et al I2OO9I : 
Ivan de Voort et al.ll2010l . hereafter VIO). 

The fraction of the accreted gas that is accreted cold, 
i.e. the fraction of the gas that did not pass through 
a virial shock, de p ends on both halo mass and redshift 
jKeres et all l2005l: lOekel fc BirnboimI I2OO6I : lOcvirk et all 
l2008l : lKeres et al.ll2009l : lBrooks et al.ll2009l : VIO). In VIO we 
used simulations to show that the rate at which gas accretes 
onto central galaxies is generally much lower than the rate at 
which gas accretes onto their host haloes. Furthermore, we 
found that while halo accretion rates are determined by the 
depth of the gravitational potentials, galaxy accretion rates 
are also sensitive to processes such as metal-line cooling and 
feedback from star formation and AGN. 

Here we use a large cosmological hydrodynamical sim- 
ulation that includes radiative cooling (computed element 
by element and thus including metal lines), star formation, 
stellar mass loss, and outflows driven by both supernovae 
and AGN, to determine how hot and cold accretion influ- 
ence the global star formation history (SFH). We calculate 
global accretion rate densities, both for hot- and cold-mode 
accretion, and for accretion onto haloes as well as accre- 
tion onto galaxies. We will compare the evolution of these 
different global accretion rates to the resulting global SFH 
and learn how they all link together. We will show that the 
sharp drop in the global SFH after z « 2 reflects the cor- 
responding sharp drop in the rate of cold-mode accretion 
onto haloes. Motivated by semi-analytic models and sim- 
ulations that have shown AGN feedback to be crucial for 
the suppression of s t ar formation in high-mass haloes (e.g . 
Sprineel et alll2005l: iBower et al.ll2006l: ICroton et al.ll200d; 



Lagos et all I2OO8I : iBooth fc Schavd l2009l : iMcCarthv et all 
20101 ). we use an identical simulation, except for the omis- 
sion of feedback from AGN, to investigate the effect of AGN 
feedback on the global accretion rates. We will show that 
the hot accretion mode is more strongly suppressed by AGN 
feedback than the cold mode. Without AGN feedback, low- 
redshift star formation would not be predominantly fuelled 
by the cold accretion mode and the drop in the cosmic SFR 
would be much less strong. 

The outline of this paper is as follows. In Section [2] we 
briefly describe the simulations from which we derive our re- 
sults and discuss our method for selecting recently accreted 
gas. In Section [3] we compare global accretion rates to the 
cosmic SFR and show which haloes contribute most to the 
global accretion rates and the cosmic SFH. In Section |3] we 
investigate the eflect of AGN feedback on the hot and cold 
accretion rate densities. Finally, we summarize and discuss 
our results in Section [S] 



2 SIMULATIONS 

We use a modifie d version of GADGET-3 (last described 
m ' Springell |2005| ). a smoothed particle hydrodynamics 
(SPH) code that uses the entropy formulation of SFH 
i Springel fc Hernauistll2002l ). which conserves both energy 
and entropy where appropriate. This work is part of 
the Overwhelming ly Large Simulations (OWLS) project 
(jSchave et al.|[2O10h , which consists of a large number of cos- 
mological simulations with varying (subgrid) physics. Here 
we make use of the so-called 'AGN' and 'reference' models, 
which are identical except that only model AGN includes su- 
permassive black holes and AGN feedback. The AGN sim- 
ulation will be our fiducial model, but we will compare it 
with the OWLS reference mo del in Section |4]_As the sim- 
ulations are fully described in ISchave et all ijioiol ). we will 
only summarize their main properties here. 

The simulations described here assume a ACDM cos- 
mology with parameters Qui = I — Qa ~ 0.238, Qh ~ 0.0418, 
h = 0.73, as = 0.74, n = 0.951 . These values are con sistenl[j 
with the WMAP year 7 data (|Komatsu et al.ll2011^ . 

A cubic volume with periodic boundary conditions is 
defined, within which the mass is distributed over iV=^ dark 
matter and as many gas particles. The box size (i.e. the 
length of a side of the simulation volume) of the simulations 
used in this work is 100 comoving Mpc, with A'^ = 512. 
The (initial) particle masses for baryons and dark matter 
are 1.2 x 10* M© and 5.6 x 10** Mq, respectively. The gravi- 
tational softening length is 7.8 /i~^comoving kpc, i.e. 1/25 of 
the mean dark matter particle separation, but we imposed a 
maximum of 2 /i~^proper kpc, which is reached sd, z — 2.91. 

The abundances of eleven elements released by mas- 
sive stars and intermediat e mass stars are followed as de- 
scribed in |w lersma et al.l (l2009bh. We assume the stellar 
initial mass function (IMF) o flChabriei'('2003^, ranging from 
0.1 to 100 Mq. As described in^iersma et al. (2009a), ra- 
diative cooling and heating are computed element by ele- 
ment in the prese nce of the cosmic microw ave background 
radiation and the iHaardt fc Madaul l|200ll ) model for the 
UV/X-ray background from galaxies and quasars. 

Star formation is mo d elled according to the recipe 
of ISchave fc Dalla Vecchial (|2008l ). A polytropic equation 
of state Ptot oc pgis is imposed for densities exceeding 
TiH = 0.1 cm~^, where Ptot is the total pressure and Pgas 
the density of the gas. Gas particles with proper densities 
riH ^ 0.1 cm"'' and temperatures T ^ 10^ K are moved onto 
this equation of state and can be converted into star parti- 
cles. The star formation rate (SFR) per unit mass depends 
on the gas pressure an d is set to repro duce the observed 
Kennicutt-Schmidt law (|Kennicuttll 19981 ) . 

Feedback from star formation is impleme nted using the 
prescription of lDalla Vecchia fc Schavd l|2008l ). About 40 per 
cent of the energy released by type II SNe is injected locally 
in kinetic form. The rest of the energy is assumed to be lost 
radiatively. The initial wind velocity is 600 km s~^. 

Our fiducial simulation includes AGN feedback. The 
mo del, which is a modifi ed version of the model introduced 
by |Springel et sT. (2005'), is fully described and tested in 
iBooth fc Schave (2009 ). In short, a seed mass black hole is 
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placed in every resolved halo. These black holes grow by ac- 
cretion of gas, which results in the injection of energy in the 
surrounding medium, and by mergers. 

The accretion rate onto the b lack hole equals the so- 
called Bondi-Hoyle accretion rate (|Bondi fc Hovlelll944l ) if 
the gas density is low (nn < 10~^ cm~^). However, in dense, 
star-forming gas, where the accretion rate would be severely 
underestimated because the simulations do not include a 
cold, interstellar gas phase and because the Jeans scales are 
unresolved, the accretion rate is multiplied by an efficiency 
parameter, a, given by a = (nu /n^i)^ , where — 0.1 cm~'^ 
is the threshold density for star formation and 13 — 2. Note, 
however, that our results are insensitive to the choice for 
P as long as it is chosen to be large (see iBooth fc Schavj 
I2OO9I ). A fraction of 1.5 per cent of the rest-mass energy of 
the accreted gas is injected into the surrounding medium in 
the form of heat, by increasing the temperature of at least 
one neighbouring gas particle by at least 10* K. The mini- 
mum temperature increase ensures that the feedback is ef- 
fective, because the radiative cooling time of the heated gas 
is sufficiently long, and results in fast outflows. When AGN 
feedback is includ ed, the SFR is reduce d for haloes with 
Mhaio > 10" M0 (|Booth fc Schave'200^. The AGN simu- 
lation reproduces the observed mass densit y in black holes at 
z = and the black hole scaling rela tions jBooth fc Schav3 
I2OO9I ) and their evolution ((Booth fc Schavd I2OI0I ') as well 
as the observed optical and X-ray properties, stellar mass 
fractions, SFRs, stellar age distribut ions and the the rmody- 
namic profiles of groups of galaxies l|McCarthv et al. 2010). 

The Lagrangian nature of the simulation is exploited by 
tracing each fluid element back in time, which is particularly 
convenient for this work, as it allows us to study the tem- 
perature history of accreted gas. During the simulations the 
maximum past temperature, Tmax, was stored in a separate 
variable. The variable was updated for each SPH particle at 
every time step for which the temperature exceeded the pre- 
vious maximum temperature. The artiflcial temperature the 
particles obtain when they are on the equation of state (i.e. 
when they are part of the unresolved multiphase ISM) was 
ignored in this process. This may cause us to underestimate 
the maximum past temperature of gas that experienced an 
accretion shock at higher densities. Ignoring such shocks is, 
however, consistent with our aims, as we are interested in 
the maximum temperature reached before the gas entered 
the galaxy. 

Resolution tests are not included in this paper. How- 
ever, ISchave et al] (|2010l ) have shown that the box size and 
resolution of the reference simulation used in this paper suf- 
fice to obtain nearly converged results for the cosmic SFH at 
z < 3, which changes by much less than a factor of 2 when 
changing the resolution by a factor of 8. At higher redshifts 
the global SFR density is, however, underestimated as a re- 
sult of the finite resolution. Because AGN feedback moves 
the peak star formation activity to lower-mass haloes, the 
convergence of the AGN simulation may be slightly worse. 
VIO tested the convergence of accretion rates and the frac- 
tion of the accretion due to the hot mode as a function of 
halo mass. They found that increasing the resolution gives 
slightly higher cold accretion fractions, which would only 
strengthen the conclusions of this work. For quantities av- 
eraged over the entire simulation volume, we have to keep 
in mind that the minimum halo mass that is included de- 



pends on the resolution. The global fraction of gas accreted 
in the cold mode may therefore also depend on the resolu- 
tion, because the cold fraction increases with decreasing halo 
mass. Increasing the resolution would allow us to include 
lower-mass haloes, thus increasing the global cold accretion 
fraction. Again, this would only strengthen our conclusions. 

2.1 Accretion and mergers 

In this section we summarize our method for determining the 
gas accretion rates onto haloes and galaxies and the amount 
of star formation within haloes. Our m ethod is described 
in more detail in VIO. We use SUBFIND (|Dolag et al.ll2009l ) 
to identify gravitationally bound haloes and subhaloes. In 
this work we only investigate accretion onto haloes that are 
not embedded inside larger haloes. In other words, we do 
not consider subhaloes. The simulations are saved at dis- 
crete output redshifts, called snapshots. The redshift inter- 
vals between snapshots are Az — 0.25 for < 2: ^ 4, and 
Az — 0.5 for 4 < z ^ 6. This is the time resolution for 
determining accretion rates. For all haloes, called descen- 
dants, with more than 100 dark matter particles (i.e. with 
A^haio ^ 10^"'* M0), we find the progenitor at the previous 
output redshift. We do this by finding the group containing 
most of the 25 most bound dark matter particles of the de- 
scendant. Haloes without a well-defined progenitor are dis- 
carded from our analysis. All gas entering a halo, both from 
accretion and from mergers, contributes to its growth. We 
therefore identify all the gas and recently formed star parti- 
cles that are in the descendant, but not in its progenitor. 

To participate in star formation, gas has to enter not 
only the halo, but also the ISM of the galaxy. For the same 
haloes as we have selected above, we identify all particles 
entering the ISM (i.e. moving onto the equation of state) 
of the descendant galaxy between the two consecutive out- 
put redshifts. Galaxy mergers are therefore automatically 
included. To make sure we are not looking at accretion onto 
substructures within the main halo, we only consider accre- 
tion taking place within the inner 15 per cent of the virial 
radius. This corresponds to the growth of the central galaxy. 

Finally, we identify all the stars formed in the same 
haloes and redshift bins to calculate the corresponding 
SFRs. Again, unresolved haloes and satellites are excluded. 

We decided to include gas contributed by mergers, as 
well as smoother accretion flows, because both are mech- 
anisms by which haloes and galaxies grow and both can 
provide fuel for star formation. We note, however, that the 
contribution from mergers is much smaller than that from 
smooth accretion. Mergers with mass ratios greater than 
1:10 account for only < 10 per cent of the total halo growth, 
except in groups and clusters at low redshift (see VIO). 

The simulations used for this study do not resolve haloes 
below Mhaio = lO^'' * Mq. We are therefore missing contri- 
butions to the global accretion rates and SFR from lower 
mass haloes. Because we use a fixed maximum past tem- 
perature threshold of 10^'^ K to distinguish hot from cold 
accretion (see below), small haloes with virial temperatures 
below this threshold will by definition accrete almost all gas 
in the cold mode. The mass regime that we investigate here 
is therefore the most interesting regime , where changes in 
accretion mode are expected to occur (|Dekel fc BirnboimI 
120061 ). 
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3 GLOBAL ACCRETION AND STAR 
FORMATION 

ISchave et al. I (|2O10l ) showed that without feedback from su- 
pernovae and AGN, the SFR density is overpredicted by 
a large factor at z < 2. By including supernova feedback 
they could lower the SFR density, but if the predicted SFR 
density matched the observed peak at z ~ 2, then it over- 
predicted the SFR density at z = 0. The drop in the global 
SFR below 2 = 2 is much closer to the observed slope if 
AGN feedback is include d, but in that case t he SFR density 
may be slightly too low (|Schave et al.ll2010l l. This discrep- 
ancy could be solved by decreasing the efficiency of feed- 
back from star formation, which had been set to reproduce 
the observed peak in the SFH using models without AGN 
feedback. It is, however, not clear how seriously we should 
take the discrepancy, since the AGN simulation does repro- 
duce the o bserved masses and a ges of the central galaxies 
of groups (jMcCarthv et al.ll2010l ). This would suggest that 
the problem may be solved by increasing as from 0.74 to 
the currently favoured value of 0.81, which would have a 
relatively strong e ffect on global star formation rates (see 
ISchave et ahlboiol ) while leaving the evolution of haloes of 
a given mass nearly unchanged. Because, for the purposes 
of this work, the AGN simulation is the most realistic run 
from the OWLS suite, we will adopt it as our fiducial model. 
In Section U we will discuss how its predictions differ from 
those of a simulation without AGN feedback. 

The temperature of accreting gas has been fo und to 
follow a bimodal distribution (e.g. iKeres et al. I l2005l : VIO). 
Therefore, two modes of accretion are considered: the hot 
and the cold mode . The cooling function peaks at about 
T = lif-^-^ K fe.g. IWiersma et al.l[2009al ). so there is little 
gas around this temperature. We therefore define hot (cold) 
accretion as accreted gas with maximum past temperatures 
above (below) 10^'^ K. Much of the gas that is shock-heated 
to much higher temperatures is expected to stay hot for a 
long time, whilst gas that is heated to lower temperatures 
can more easily cool and condense onto the galaxy. 

When considering accretion onto galaxies, it is impor- 
tant to note that the terms 'hot' and 'cold' refer to the maxi- 
mum past temperature rather than the current temperature. 
In fact, gas that has been shock-heated to temperatures in 
excess of 10^'^ K must cool down before it can accrete onto a 
galaxy, as we identify this type of accretion as gas joining the 
ISM, which in our simulations means that the gas density 
crosses the threshold of nn =0.1 cm~^ while the tempera- 
ture T ^ 10^ K. Thus, although gas that, according to our 
terminology (which is consistent with that used in the litera- 
ture), is accreted hot onto galaxies has been hot in the past, 
it is cold at the time of accretion. In VIO we showed that 
the maximum past temperature is usually reached around 
the time the gas accreted onto the halo. 

The global accretion rate density onto haloes, i.e. the 
gas mass accreting onto resolved haloes per year and per 
comoving Mpc^, is shown in the left panel of Figure [T] by 
the solid, black curve. The solid, red and blue curves show 
global accretion rates for hot- and cold-mode accretion, re- 
spectively. These global accretion rate densities are averaged 
over the time interval between two snapshots. The SFR den- 
sity in the haloes we consider is shown by the green curve. 
The SFH in the entire box is shown by the dashed, grey 



curve. It is higher than the dashed, green curve, because the 
latter excludes star formation in subhaloes and unresolved 
haloes (i.e. Mhaio < lO^"'* Mq). 

The global accretion rate onto resolved main haloes 
(solid, black curve) peaks at z f» 3. It is fairly constant, 
varying only by about a factor of two from z « 4 down to 
z ~ 0. The average accretion rate onto haloes of a given 
mass decreases more strongly towards lower redshift (VIO). 
However, the number of haloes at a fixed mass increases and 
higher mass haloes form with decreasing redshift. The com- 
bination of these effects results in an almost constant global 
accretion rate density. We note, however, that the normal- 
ization and shape are not fully converged with respect to res- 
olution: higher resolution simulations may find higher total 
and cold accretion rate densities. Increasing the resolution 
would allow us to include haloes with lower masses, which 
would boost the cold accretion rate density by up to a factor 
of ~ 2 — 3. The total accretion rate density would there- 
fore increase at z > 2, where cold accretion dominates. The 
global accretion rate is an order of magnitude higher than 
the global SFR in the same haloes (green, dashed curve), in- 
dicating that most of the gas that accretes onto haloes never 
forms stars. 

The global growth of haloes is dominated by cold ac- 
cretion for z > 2, but by z = the contribution of the hot 
mode exceeds that of the cold mode by an order of magni- 
tude. The cold accretion rate density peaks around z = 3 
and falls rapidly thereafter, while the hot accretion rate den- 
sity increases down to z ~ 2 and fiattens off at lower red- 
shifts. The global SFR peaks around z = 2 and declines by 
an order of magnitude to z = 0. An order of magnitude drop 
is also visible for the global cold accretion rate from z — 3 
to z = 0. Neither the total nor the hot accretion rate his- 
tories can explain the drop in the cosmic SFH, which must 
therefore be driven by the drop in the cold accretion rate. 

The middle panel of Figure [T] shows the global accretion 
rates onto central galaxies. The green and grey curves are 
identical to the ones shown in the left panel. The black, red, 
and blue curves describe the total, hot, and cold accretion 
rate densities, respectively. The total accretion rate shows 
the amount of gas that joins the ISM. The gas can, however, 
be removed from the ISM by supernova and AGN feedback, 
as well as by dynamical processes. This is why the overall 
normalization of the SFR is generally lower than the ISM 
accretion rate. The global SFR peaks later than the global 
galaxy accretion rate (z « 2 versus z « 3, which corresponds 
to a difference of ~ 1 Gyr in the age of the Universe). This 
delay probably results from the time it takes to convert 
interstellar gas into stars. The long gas consumption time 
scale implied by the assumed, Kennicutt star formation law 
(~ 1 Gyr for typical densities) allows the existence of reser- 
voirs of accreted gas. The SFR can therefore temporarily be 
higher than the galaxy accretion rate, as happens for z < 0.5 
(compare the dashed, green curve with the solid, black curve, 
which include the same sample of resolved haloes). Gas re- 
turned to the ISM by stellar mass loss, a process that is 
include d in our simulations and that becomes importa nt for 
z < 1 (jSchave et al.ll2010l : [Leitner fc Kravtsovll2010l l. also 
increases the SFR relative to the accretion rate. 

As we did in the left panel of Figure [T] for accretion onto 
haloes, we split the global accretion rate onto galaxies into 
separate contributions from the hot and cold modes. The 
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Figure 1. Evolution of global accretion rate and SFR densities in resolved haloes with well-defined progenitors. The left, middle, and 
right panels show global accretion rate densities onto haloes, onto central galaxies, and global SFR densities, respectively. Red and blue 
curves show accretion rate densities (left and middle panels) and SFR densities (right panel) resulting from hot- and cold-mode accretion, 
respectively. In all panels, the black curve is the sum of the red and blue curves, the green curve shows the global SFR density in the 
selected haloes, and the grey curve shows the global SFR density in the entire simulation box. The small 'step' visible at 2 ft: 4 is caused 
by the sudden increase in the time resolution for determining accretion, i.e. Az between snapshots decreases by a factor of two at 2 = 4. 
Galaxies accrete most of their gas in the cold mode and this mode is responsible for an even larger fraction of the star formation. Because 
of outflows driven by supernovae and AGN, the SFR density is generally lower than the galaxy accretion rate density. The global SFR 
declines more rapidly than either the total and hot-mode accretion rate densities. This decline must therefore be caused by the drop in 
global cold-mode accretion rate, though with a delay. 



global hot accretion rate peaks around z = 2, as does the 
SFR density. The hot accretion rate is, however, not nearly 
enough to account for all the star formation in these galaxies, 
falling short by at least a factor of 2 at all redshifts. At z > 3, 
the global cold accretion rate is an order of magnitude higher 
than the global hot accretion rate. This difference decreases 
to ~ 0.25 dex by z = 0. At all redshifts, it is mostly cold 
accretion that allows for the growth of galaxies, even though 
hot accretion dominates the growth of haloes below 2: ~ 2. 

Comparing the middle panel to the left panel, we notice 
that the global cold accretion rate onto galaxies is a factor of 
~ 3—4 lower than the cold accretion rate onto haloes. Not all 
cold gas that accretes onto haloes makes it into the central 
galaxy to form stars (see also VIO). The shapes of the blue 
curves are, however, similar, indicating that the fraction of 
the gas accreting cold onto the halo that proceeds to accrete 
onto the galaxy is roughly constant with time. 

The situation is very different for global hot accretion 
rates (red curves). Hot accretion onto the ISM has already 
peaked at z « 2.5, while hot accretion onto haloes con- 
tinues to increase down to z = 0. This can be explained 
by noting that as the Universe evolves, gas is heated to 
higher temperatures, because haloes become more massive. 
In addition, the average density of the Universe goes down. 
The lower densities and higher temperatures gives rise to 
longer cooling times. Moreover, winds from supernovae and 
AGN eject low-entropy gas at high redsh ift, raising the en- 
tropy of the gas in ha loes at low redshift (ICrain et al.ll201ol : 
iMcCarthv et al.ll201ll ). Hence, as time goes on, more of the 
halo gas is unable to cool and reach the central galaxy. While 
the gravitational potential is the most important factor for 
the growth of haloes, for the growth of galaxies, the cooling 
function and feedback processes also come into play. 

The right panel of Figure [U which shows the global 
SFR densities due to hot and cold accretion, confirms that 
the main fuel for star formation is gas accreted in the cold 



mode. The difference becomes smaller towards lower red- 
shift. However, even at z = hot mode gas contributes 
0.3 dex less than cold mode gas. 

To investigate which haloes contribute most to the 
global accretion rates and SFHs, we show the same quan- 
tities as in Figure [1] for three different halo mass bins 
in Figure [2] From top to bottom, the mass ranges are 
10" < Mhaio < 10^2 Mq, 10^2 < Mhaio < 10'=^ M©, and 
Mhaio ^ 10^^ Mq, which contain 21813, 2804, and 285 haloes 
at z = 0. The shape of the total halo accretion rate density 
in the lowest mass r ange is in agreement with that found by 
iBouche et all (j20ld ) based on an extended Press-Schechter 
formalism and a fit to dark matter accretion rates in N-body 
simulations. 

For z > 2 the global accretion rate densities onto haloes 
(left column), galaxies (middle column), and the global star 
formation rate density (right column) are all dominated by 
haloes with Mhaio < lO^'^ Mq (top row). At that time, 
higher-mass haloes are still too rare to contribute signif- 
icantly. Below z « 2 haloes with Mhaio = 10^^"^^ Mq 
(middle row) begin to contribute significantly and for ac- 
cretion onto haloes, but not for accretion onto galaxies or 
star formation, their contribution is overtaken by Mhaio = 
]^Qi3-i4 haloes (bottom row) around 2 = 1. 

Observe that the global halo accretion rate density 
starts to decline at z « 3 and z « 2 for the low and mid- 
dle halo mass bins, respectively, and that it keeps increasing 
down to z = for Mhaio > 10^^ Mq. The global cold ac- 
cretion rate density decreases with time for z < 3, z < 2.5, 
and z < 1 for the low, middle, and high halo mass bins, 
respectively. 

Both the galaxy accretion rate density and the SFR 
density are dominated by Afhaio < lO'^^ Mq haloes at all 
redshifts. Towards lower redshifts, high-mass haloes account 
for larger fractions of the total galaxy accretion rate density 
and SFR density, though they never dominate. High-mass 
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Figure 2. Evolution of global accretion rate densities onto haloes (left column), galaxies (middle column), and global SFR densities 
(right column) for different halo mass bins. From top to bottom only haloes have been included with: 10^ ^ My^^i^ < 10^^ M©, 
10^^ Afhalo < lO'^^ Mq, and Mij^io ^ lO'^^ Mq. The curves show the same quantities as in Figure[T] Above z fa 3.5 the highest mass 
bin contains no haloes. At all redshifts most of the cold halo accretion, galaxy accretion, and star formation happens in low-mass haloes 
(i.e. Mhaio < 10^2 Mq). At z > 2 low -mass haloes also dominate the global hot halo accretion rates. At 2 < 1 the total halo accretion 
rate is dominated by high-mass haloes, but nearly all of the gas is accreted hot and unable to accrete onto galaxies. 



haloes do dominate the halo accretion rate at low redshift, 
but nearly all of the gas is accreted hot and only a very small 
fraction of this gas is subsequently able to cool down onto 
galaxies. 



4 EFFECT OF AGN FEEDBACK 

It is interesting to see what the influence of AGN feedback 
is on our results. Because AGN feedback is more important 
in higher-mass haloes, for which hot-mode accretion is more 
important, we expect it to have a larger effect on the hot- 
mode accre tion rate density. Moreover, it has been hy poth- 
esized (e.g. iKeres et al.]|2005l : iDekel fc Birnboimllioog ') that 
hot, dilute gas may be more vulnerable to AGN feedback 
than cold streams and may therefore be pref erentially pre- 
vented from accreting onto galaxies. Indeed, Ixheuns et al.l 
l|2002l) had already demonstrated that supernova-driven out- 



flows follow the path of least resistance, leaving the cold 
filaments that produce HI absorption intact. 

iMcCarthv et al.l (|201ll ) have shown that feedback from 
AGN at high redshift increases the entropy of the halo gas 
at low redshift. The hot gas will therefore be even hotter and 
less dense at low redshift than it would be in the absence of 
AGN feedback, making it more susceptible to being heated 
or entrained in an outfiow, and thus to being prevented from 
accreting. 

We compare our fiducial simulation, which includes 
AGN feedback, to the OWLS 'reference model' which is 
identical to our fiducial run except that it does not include 
black holes and AGN feedback. This allows us to assess the 
effect of AGN feedback on the global hot and cold accretion 
rates. Figure shows the same solid curves as were shown 
in Figure [1] They indicate the total, hot, and cold accre- 
tion rate densities onto haloes (left panel) and onto galaxies 
(middle panel) and the star formation rate density resulting 
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Figure 3. Evolution of global accretion rate densities onto haloes (left column), galaxies (middle column), and global SFR densities for 
simulations with (solid curves, same as in Figure[lJ and without (dashed curves) AGN feedback. The black, red, and blue curves show the 
global accretion and SFR densities from all, hot-mode, and cold-mode accretion, respectively. AGN feedback suppresses halo accretion 
only slightly, but the effect on galaxy accretion and star formation is large, up to an order of magnitude. AGN feedback preferentially 
suppresses hot-mode galaxy accretion and star formation from gas accreted in the hot mode. 



from all, hot, and cold accretion (right panel). The dashed 
curves show the same global accretion rates and SFHs for 
the simulation without AGN feedback. For accretion onto 
galaxies and for star formation the differences are striking. 
When AGN feedback is excluded, late-time star formation is 
no longer predominantly fuelled by gas accreted in the cold 
mode. 

As expected, all accretion rate densities are reduced by 
the inclusion of AGN feedback. The effect on halo accretion 
is, however, small, as was also shown by VIO. The hot and 
cold halo accretion rate densities are reduced by at most 
0.2 and 0.1 dex, respectively. This reduction implies that 
AGN feedback also affects some gas outside of haloes. Even 
though the effect is small, AGN feedback reduces hot halo 
accretion more than cold halo accretion. 

The differential effect of AGN feedback on hot and cold 
accretion is much more pronounced for accretion onto galax- 
ies than for accretion onto haloes and it incre ases towards 
lower redshift. At very high redshift {z = 9). I Powell et al.l 
(|2010l ) have shown that outflows (driven by supernova feed- 
back) do not affect the galaxy inflow rates. Our results in- 
dicate that this may change towards lower redshifts, when 
densities are much lower. While AGN feedback reduces cold 
accretion rate densities by up to 0.4 dex (at z — 0), the hot 
accretion rate densities decrease by up to 0.8 dex (also at 
z — 0). The SFR densities are reduced by up to 0.6 dex for 
star formation powered by cold-mode accretion, but by 1 dex 
for hot-mode accretion. The reduction due to AGN feedback 
is thus ~ 0.4 dex greater for the hot mode than for the cold 
mode, both for galaxy accretion and for star formation. The 
larger reduction indicates that AGN feedback preferentially, 
but not exclusively, prevents hot mode gas from accreting 
onto galaxies and participating in star formation. 

Hence, the inclusion of AGN feedback strongly boosts 
the size of the drop in the cosmic SFR at late times. This 
preferential suppression of hot accretion is the result of two 
effects, namely of the differential effect at a flxed halo mass, 
indicating that hot-mode gas is more vulnerable to feedback 
than cold-mode gas, and of the fact that AGN feedback is 
effective only in massive haloes (with Afhaio 10^^), for 



which hot accretion is important. The latter is the dominant 
effect. 



5 CONCLUSIONS 

We have investigated the evolution of the global gas accre- 
tion rate densities onto haloes and onto their central galax- 
ies and we have done so for both the hot and cold accretion 
modes. In addition, we studied the contributions from gas 
accreted through the cold and hot modes to the cosmic star 
formation history. We made use of a 100 Mpc/h, 2x512'^ par- 
ticle SPH simulation from the OWLS project that includes 
radiative cooling (computed element by element and thus 
including metal lines), star formation, stellar mass loss, su- 
pernova feedback, and AGN feedback. We isolated the effect 
of AGN feedback by comparing to a second simulation that 
did not include AGN, but which was otherwise identical. The 
hot and cold accretion modes were separated by using a flxed 
maximum past temperature threshold of Tmax ~ 10^'^ K. 

The global gas accretion rate density onto haloes is 
much higher than that onto galaxies and both rates exceed 
the cosmic SFR density. This confirms the finding of VIO 
that most of the gas accreting onto haloes does not result 
in star formation. This is the case for both accretion modes, 
but the differences are larger for the hot mode. 

The global SFR declines after z ~ 2, whereas the global 
hot-mode accretion rate onto haloes shows no such trend. 
From this, we conclude that the global SFR follows the drop 
in the global cold-mode accretion rate onto haloes, which 
sets in at z r; 3, but with a delay of order the gas con- 
sumption time scale in the ISM. Star formation tracks cold- 
mode accretion rather than hot-mode accretion because cold 
streams can reach the central galaxy, where star formation 
takes place, much more easily than gas that is shock-heated 
to high temperatures near the virial radius. Much of the hot 
gas cannot cool within a Hubble time and therefore can- 
not accrete onto the central galaxy. In addition, we demon- 
strated that it is very important that hot gas is more suscep- 
tible to removal by outflows driven by feedback from AGN. 
Without AGN feedback, gas accreted in the hot mode con- 
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tributes significantly to the cosmic SFR below z = 1 and 
the drop in the SFR below z = 2 would be much smaller. 

For the hot mode the difference between the accretion 
rates onto haloes and onto galaxies is larger at lower red- 
shifts. While the hot accretion mode dominates the growth 
of haloes by an order of magnitude at z « 0, it is still less 
important than cold accretion for the growth of the central 
galaxies. At z > 2, cold accretion even dominates the global 
accretion rate onto haloes. 

We demonstrated that AGN feedback suppresses accre- 
tion onto galaxies and that it does so much more efficiently 
for the hot mode than for the cold mode. This happens 
because AGN feedback only becomes more efflcient than 
feedback from star formation in high-mass haloes, which 
are dominated by hot accretion, and because hot-mode gas 
is more dilute and therefore m ore vulnerable to feedback. 
In addition, as demonstrated bv lMcCarthv et al.l (|201ll l. by 
ejecting low-entropy halo gas at high redshift (z > 2), AGN 
feedback results in an increase of the entropy, and thus a re- 
duction of the cooling rates, o f hot halo gas at low redshift. 

While iKeres et al.l (|2009l ) did not investigate accretion 
onto haloes, they did also find that cold accretion is most im- 
portant for the growth of galaxies, with hot accretion becom- 
ing increasingly important towards lower redshifts (see also 
iKeres et al.l l2005l : locvirk et all l2008l : iBrooks et all l2009l : 
VIO). However, their simulation included neither winds from 
supernovae nor feedback from AGN. AGN feedback was in 
fact ignored by all previous cosmologi cal simula tions investi- 
gating gas accretion except for Khal atvan et al., (2008), who 
simulated a single object, and except for VIO. Our results 
suggest that the neglect of this important process leads to 
a strong overestimate of the global accretion rate and SFR 
densities and of the importance of the hot accretion mode 
for galaxy accretion and star formation. 

In summary, the rapid decline in the cosmic SFR density 
below z = 2 is driven by the corresponding drop in the cold 
accretion rate density onto haloes. The total accretion rate 
onto haloes falls off much less rapidly because the hot mode 
becomes increasingly important. AGN feedback, which acts 
preferentially on gas accreted in the hot mode, prevents the 
hot halo gas from accreting onto galaxies and forming stars 
and is therefore a crucial factor in the steep decline of the 
cosmic SFR density. 
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